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Abstract: 

We present an algorithm to compute arbitrary multi-loop massive Feynman diagrams 
in the region where the typical energy scale y/s is much larger than the typical mass scale 
M, i.e. s ^ M^, while various different energy and mass parameters may be present. 
In this region we perform an asymptotic expansion and, using sector decomposition, we 
extract the leading contributions resulting from ultraviolet and mass singularities, which 
consist of large logarithms ln(s/M^) and 1/e poles in D = 4 — 2e dimensions. To next-to- 
leading accuracy, at L loops all terms of the form a^e~'^ In-' (s/M^) with j + k = 2L and 
j + k = 2L — 1 are taken into account. This algorithm permits, in particular, to compute 
higher-order next-to-leading logarithmic electroweak corrections for processes involving 
various kinematical invariants of the order of hundreds of GeV and masses Mw ~ Mz ~ 
Mh ~ mt of the order of the electroweak scale, in the approximation where the masses of 
the light fermions are neglected. 
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1 Introduction 



With the advent of colliders in the TeV energy range, multi-loop calculations, both in 
QCD and in the electroweak theory, will play an increasingly important role in order to 
provide theoretical predictions of sufficiently high accuracy. On the one hand, this is due 
to the high level of experimental precision, especially at the linear collider [1-3]. On the 
other hand, the size of radiative corrections grows with the energy. 

This latter feature is particularly important for the electroweak corrections [4], since 
at TeV colliders the energy scale ^/s starts to be much larger than the weak-boson mass 
scale M = Mw ^ Mz- For <C s, i.e. if one performs an asymptotic expansion in 
M^/s, the electroweak corrections assume the form of a tower of large logarithms. 



at L-loop level. These logarithms originate from ultraviolet (UV) and mass singularities. 
Above the electroweak scale, such logarithms start to grow and become the leading con- 
tribution to the electroweak corrections. At ^/s = 1 TeV they give rise to corrections of 
tens of per cent at one loop and a few per cent at two loops. Thus, in order to reach the 
per-mille level accuracy envisaged by the linear collider, it becomes mandatory to compute 
electroweak radiative corrections up to the two-loop level. 

Despite of recent progress in this direction [5], the exact analytical computation of 
two-loop electroweak corrections for processes involving many (more than three) external 
legs and various internal masses is out of sight. One is therefore forced either to adopt a 
numerical approach, as for instance in Refs. [6,7], or to consider an asymptotic expansion 
of the type (1), as we do in the present paper. 

The asymptotic behaviour of higher-order electroweak corrections has received a lot of 
interest in the recent years [8-21]. On the one hand, resummation prescriptions have been 
proposed [8-12], which predict the leading logarithms (LLs) and next-to-leading logarithms 
(NLLs), i.e. terms with j = 2L and 2L — 1 in (1), for arbitrary processes and also the 
next-to-next-to-leading logarithms (NNLLs), with j = 2L — 2, for 2 — >^ 2 massless fermionic 
processes [13]. On the other hand, in order to check these resummations, which rely on the 
assumption that various aspects of the symmetry-breaking mechanism can be neglected in 
the high-energy limit, explicit diagrammatic two-loop calculations have been performed. 
For the LLs [15-17] as well as for the angular-dependent subset of the NLLs [18] such 
checks have been already completed and the exponentiation predicted in Refs. [10 -12] has 
been confirmed for arbitrary processes. At (or beyond) the level of the NLLs only very 
few calculations exist. In the electroweak Standard Model, a calculation of a massless 
fermionic form factor [19] confirmed the resummation prescriptions of Refs. [11,12] up to 
the NLLs. The evaluation of the complete tower of logarithms contributing to this form 
factor in the special case of a U(l) x U(l) theory with a massive and a massless gauge 
boson [20,21] was found to be in agreement with the resummation prescriptions of Ref. [13] 
at the NNLL level. 

For general processes, the terms withj = 2L— 1 in the asymptotic expansion (1) remain 

largely untested and those for j < 2L — 2 almost completely unknown. It is important to 
further investigate the behaviour of these subleading logarithms since, depending on the 
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process, their phenomcnological impact at the TeV scale might be comparable with that 
of the LLs, as was found in Refs. [11,13,20,21] for processes involving massless fermions. 

In this paper we present an algorithm to compute the LLs and NLLs in a fully au- 
tomated way. This algorithm, which has already been used in Refs. [18,19], applies to 
arbitrary multi-loop Fcynman diagrams involving various energy parameters Sj ~ s and 
masses ^ M oi the same order. In the asymptotic region s 3> M^, the coefficients of 
the LLs and NLLs are computed analytically as a function of the ratios sj/s and M^/M. 
The mass singularities originating from soft/collinear massless particles (photons, gluons, 
or fermions) as well as the UV singularities are regulated dimensionally and appear as 1/e 
poles in D = 4 — 2£ dimensions. 

The algorithm is based on the so-called sector decomposition, a technique that permits 
to separate overlapping UV or mass singularities in Feynman-parameter (FP) integrals. 
Originally, this method was introduced in order to deal with UV singularities [22]. Later, 
it turned out to be applicable to mass singularities [23,24], and a first general algorithm to 
isolate 1/e poles in massless loop integrals has been formulated in Ref. [7] . The algorithm 
that we present here permits to extract combinations (products) of 1/e poles and mass- 
singular logarithms of s/M^, which arise in the case where both massless and massive 
particles are present. 

The input is an arbitrary multi-loop tensor integral, which is first translated into a 
FP integral. Depending on the location of UV and mass singularities in FP space, the 
integration range is decomposed into sectors, and each sector is remapped into the original 
integration range. This can be done iteratively until, in each sector, all singularities arise 
when individual FPs tend to zero and these FPs can be factorized. As a consequence, 
a simple power counting permits to determine the degree of singularity, i.e. the power of 
logarithms and l/e poles resulting from each sector. One can thus easily select those sectors 
that are responsible for the leading and next-to-leading singularities. Furthermore, those 
FP integrations that lead to singularities have a standard structure and can be computed 
in analytic form. 

The coefficients of the resulting logarithms and l/e poles are represented as integrals 
over the remaining FPs. As a result of sector decomposition, these last integrations are 
free of singularities and, within the Euclidean region, the integrands are smooth functions 
which can be integrated numerically. Moreover, the coefficients of the leading and next-to- 
leading singularities arc very simple integrals and in all examples we considered up to now 
(two-loop diagrams with up to 4 external legs and various internal masses) they could be 
solved in analytic form by standard computer algebra programs such as Mathematica. 

The paper is organized as follows: in Sect. 2 we describe the basics of the high-energy 
limit and the logarithmic approximation. In Sect. 3 we review the general expressions for 
multi-loop FP integrals and their UV and mass singularities and sketch the idea of the 
sector decomposition. The sector decomposition and the extraction of the singularities 
of massless diagrams is presented in Sect. 4 and the treatment of massive diagrams in 
Sect. 5. The applicability of the proposed method is discussed in Sect. 6. Some notations 
and conventions and several results for relevant massless and massive integrals are listed 
in the appendices. 
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2 High-energy limit and logarithmic approximation 



The computation of higher-order electroweak corrections is a multi-scale problem. Typi- 
cally there are various energy parameters si,S2, ■ ■ ■, which correspond to the kinematical 
invariants of the process, masses of the order of the electroweak scale M = Mw ~ Mz ~ 
Mh ~ rrit, various light-fermion masses <C M and a fictitious infinitesimal photon 
mass to regulate infrared singularities. If all light fermions (including b quarks) and pho- 
tons are treated as massless, then all particle masses arc cither zero or of the order M. 
Assuming this approach, in the present paper we consider the kinematical regime where all 
energy parameters are much higher than the electroweak scale. To be general, we consider 
multi-loop Feynman diagrams characterized by the hierarchy 

s ~ |si| ~ |s2| ~ . . . ~ \sj\ M"^ ^ M2 ^ . . . ^ Ml; (2) 

of energy parameters and masses, i.e. we assume that there are only two different scales 
s > and > 0. The external particles can be either on-shell with masses of the order 
M or zero, or off-shell with invariant masses of the order s. In this regime, we perform an 
asymptotic expansion in the small mass-to-energy ratio 

= — < 1, (3) 

treating the ratios Sj/s and M^/M as constants in the limit w ^ 0. The leading terms of 
this expansion arc divergent. On the one hand there arc mass singularities associated to 
massive particles which appear as logarithms of w. On the other hand, mass singularities 
from massless particles and UV singularities appear as 1/e poles in D = 4 — 2e dimensions. 
The leading terms of the asymptotic expansion of L-loop Feynman integrals^ assume the 
form of a double series in e and ln(u;), 

2L 00 

A{w) =w'^J2Yl "j-.fe ^'^ ^^^^Hw) + 0{w'^+^). (4) 

j=0 k=-j 

In general, the loop integrals can give rise to non-logarithmic mass singularities, i.e. d can 
be negative. However, in spontaneously broken gauge theories such singularities are com- 
pensated by mass terms or coupling constants proportional to masses in such a way that 
the resulting contributions are at most logarithmically divergent. In the present paper, 
we compute loop integrals restricting ourselves to the leading terms of order w'^ whereas 
contributions of order w'^'^^, which are suppressed by powers of masses, are neglected. 

The various terms in (4) are classified according to their degree of singularity j, which 
corresponds to the total power of logarithms and 1/e poles. For a fixed value of j we 
compute all terms of the order 

e-^,e-^+^ ln{w), 6'^+^ \n'^{w), {w),e\n^+\w), (5) 

up to the needed order in e. At L-loop level, all terms of the type (5) with j = 
2L, 2L — 1, 2L — 2, . . ., are collectively denoted as leading logarithms (LLs), next-to-leading 

^We assume that overall factors with mass dimension have been factorized such that the ampUtude 
A{w) is dimensionless. 
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logarithms (NLLs), next-to-ncxt-to-leading logarithms (NNLLs), and so on. For instance, 
with this convention, the logarithms of order In^^^w) as well as the poles of order e~'^^ 
are denoted as LLs. In the present paper we adopt a NLL approximation, i.e. we include 
all contribution of the type (5) with j >2L — 1. To denote this approximation we use the 

, 1 NLL 

symbol = . 

Since there are different energies and masses, the coefficients of the expansion (4) 
depend on their ratios, sj/s and M^/M. This dependence is taken into account and 
permits, for instance, to compute electroweak corrections without assuming that Mw = 
Mz = Mh = mt and/or including angular-dependent corrections. 



3 Fey nman- parameter integrals 

In this section we review the well-known technique [25] to translate loop integrals into FP 
integrals. Then we discuss the origin of UV and mass singularities in FP space. 



3.1 Derivation of the Feynman-parameter representation 

Let us consider a generic L-loop Fcynman integral in D dimensions with / internal lines, 
E external lines, and m loop momenta in the numerator, 

nd«|-^ d«-M<'-'^. (6) 

The propagators Dj^ = (r^ — m^ + ie)^"'^ with infinitesimal positive imaginary part e, which 
are raised to integer powers > 1, are determined by the internal momenta r^, which 
consist of linear combinations of loop momenta q = (gi, . . . , g^) and external momenta p = 
(pi, . . . ,pe), depending on the topology. The loop momenta occurring in the numerator are 
contracted with a projector -R^^...^^, which depends in general on the external momenta. 
In the following, we choose to split off a universal factor (47r/i^/s)*^^^^/^'''^ from the L- 
loop Feynman integral in order to simplify the notation. This is equivalent to the choice 
//^ = s/(47r). This factor can be trivially restored at the end if needed. 

In order to combine the products of propagators into a single denominator one intro- 
duces FPs by means of the well-known identity 

with = Xls=i "s- III order to perform the integration over the loop momenta q, we 
eliminate the linear terms from the quadratic form in the denominator, so that 

ILL L 

asDs = Qi^ijQl - 2 ^ qiQi - J + ie = J2 hMijh - A + ie. (8) 
This is easily obtained by a shift of the loop momenta with 

L L 
ki = qi- VI, vi = J2 M^^'Qu A = J+J2 Q^M-^Qj. (9) 
i=l i,j=l 
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The resulting FP integral consists of a sum of terms with different numbers of loop mo- 
menta in the numerator. Only terms with an even number of loop momenta do not vanish, 
i.e. terms of the type 



G 



T{N) 



nf=ir(n. 



•^0 ^ s=l ^ s=l 



X 



i-Mi . . . 



tl2p+l o 



(10) 



with < 2p < m. The integration over the loop momenta yields 

L 



r(iv) / \{dki 

1=1 



"■ai • • • '^0'2p 



[U^Mk - A + ie 



N 



IS 



2-D/2 



(47r)2 



(det M) 



-D/2 



where 7(M ^) = 1 for p = 0, whereas for p > 1 



^^ ■■■In h=A « -^f^ 



(11) 



(12) 



il...ip fc=l 



with the sum running over all (2p)!/(p! 2^) pairings (ii, ji), . . . , {ip,jp) obtained from the 
set {1,2,..., 2p}. As a result we are left with a FP integral of the form 



G = ^d'ddi^l-Y^as^m- 



In order to isolate the inverse powers of detM which originate from M ^, we write the 
integrand g{d) in terms of ^ 

W = detM, M.~^ = M.~Met M, = s~^Adet M - ie, Wi = VidetM, (14) 

resulting in 

Af{d) 



g(d) = T(e + Le) ,^ r /r,1^ 

[J^{d)]^+^' [uid)]f-^^+^^' 

in D = 4 — 2£ dimensions, with e = N — 2L — p, f = 2(L + 1) — N + m and 



(15) 



Mid) = {-i)^-Ps-^ 



(4vr) 



n 



Ta^'.Vaa^M-^) • • • i?^,...^.. (16) 



Note that we factorize the energy scale s in order to define as a dimensionless function. 
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3.2 Ultraviolet and mass singularities 

As discussed in Sect. 2, the large logarithmic contributions that dominate the high-energy 
expansion of loop diagrams originate from UV and mass singularities. Let us now discuss 
the origin of such singularities in the FP integrals (13)-(16). 

The overall gamma function r(e + e) in (15) represents an obvious source of UV singu- 
larities if e < 0. In addition, there are less trivial UV and mass singularities that originate 
from end-point singularities of the integrand. Let us focus on this kind of divergences, 
which result from the zeros of the polynomials lA and/or T when their exponents, / and 
e, are positive. 

The zeros of the function lA, which are independent of the external masses and mo- 
menta, are associated to UV sub divergences. Instead, the zeros of the function T can 
always be avoided by choosing non-zero internal and external masses and thus correspond 
to mass singularities. These latter, as is well known, originate from soft and/or collinear 
regions in momentum space. 

It is useful to recall that the functions U and which are independent of the loop 
momenta occurring in the numerator of (6), are determined by the topology of the corre- 
sponding Feynman diagram [25]. First, U is an homogeneous polynomial of degree L in 
the FPs, 

U{d) = Y^X{ai,, (17) 

T k=l 

where the sum runs over all trees T of the Feynman diagram. A tree is obtained by cutting 
L lines (zi, . . . ,zl) of the diagram such that all vertices remain connected. Second, is a 
homogeneous polynomial of degree L -|- 1 in the FPs, 

L+l / 

s J"(a) = - ^se ai^+U{a)^ cxsml - ie, (18) 

C k=l s=l 

where the first sum runs over all cuts C of the Feynman diagram. A cut is a set of L-|-l lines 
(ii, . . . ,iL+i) that when cut divide the diagram in two connected subdiagrams. Finally, 

sc = Vi^^ with fi^ = ±ri^ (19) 

denotes the squared total momentum^ flowing through the cut C. 

All coefficients of the polynomial U are non-negative. This implies that the UV sub- 
divergences, which result from the zeros of this polynomial, originate only from regions of 
the type 

{a\aj^= ... = aj„ = 0} with 1 < n < I - 1, (20) 

where one or more FPs vanish. Here {ji, . . . ,jn} is a subset of {1, . . . ,/}. Mass singular- 
ities originate from the zeros of the polynomial ^ in the limit of vanishing internal and 
external masses, where T depends only on the invariants sc that are not squares of on-shell 

'^Thc signs in (19) must be chosen in such a way that all momenta fj^ have the same direction with 
respect to the cut C. 
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external momenta, i.e. the invariants .s^-. In the case where these invariants arc negative 
{sj < 0) all coefficients of the polynomial (18) arc non- negative for vanishing internal and 
external masses, and mass singularities only originate from regions of the type (20). This 
is valid for arbitrary, i.e. also positive, values of the invariants sj since mass singularities 
appear independently of the directions of the external momenta [26] and thus of the values 
(and the signs) of these invariants. 



3.3 Subtraction of overlapping singularities 



Using the fact that all singularities we are interested in originate from regions of the type 
(20), we can isolate them at the level of the integrand and simplify their computation. 
This can be done with simple subtractions that split the integrand into singular and 
non-singular parts in such a way that the integration of the singular part is drastically 
simplified. 

Let us first illustrate how to proceed in the simple case of a one-dimensional integral 
assuming that the integrand has the form l/[a/o(a)] with a polynomial /q such that 
/o(a) > for < a < 1. Here, the logarithmic singularity can be isolated by means of a 
simple subtraction at a = 0: 



/■I da 1 _ /"^ da 1 _^ da 
Jo a /o(a) Jo a /o(0) Jo a 



1 



1 



/o(a) /o(0) 



(21) 



singular 



finite 



As a result, the singular part, which requires the introduction of an appropriate regular- 
ization prescription, becomes much simpler than the original integral. Also the remaining 
non-singular part is simplified since the integrand is a smooth function and can be inte- 
grated, either analytically or numerically, in the absence of a regulator. 

In the case of massive integrals, the mass-to-energy ratio w = M'^/s plays the role of 
the regulator and the mass-singular logarithms that arise in the high-energy limit, w <C 1, 
can be isolated by means of a subtraction similar to (21). If we assume that the integrand 
has the form l/[a/o(a) + wfi{a)], with polynomials fo and /i such that fi{a) > for 
< a < 1, then we can write 



da 



da 



Jo afo{a) +wfi{a) Jo afo{0) + wfi{0) 



JO a 



1 



1 



Ma) /o(0) 



+ 0{w), (22) 



and the mass-suppressed terms of order w can be neglected. 

In order to be able to perform the above subtractions, it is crucial that the polynomial 
which leads to the singularity is written in the form a/o(a), where the integration variable 
a, which is responsible for the singularity at a = 0, is factorized. In the above one- 
dimensional examples this is always possible, since any polynomial with a zero at a = can 
be written in this form. Instead, in the case of multi-dimensional FP integrals, the situation 
is complicated by the presence of so-called overlapping singularities, which originate only 
from regions where different FPs vanish simultaneously, i.e. regions of the type (20) with 
n > 1. The polynomials that are responsible for these overlapping singularities can be 
written in the form J2k=i aj^.fo,kia) and the FPs a^j. cannot be factorized. 

This difficulty can be overcome by means of the so-called sector-decomposition tech- 
nique, which consists of a decomposition of the integration range into various sectors 
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followed by a remapping of each sector into the original integration range. As we show 
in the following two sections, the sectors can be chosen in such a way that the resulting 
integrals are of the form 

(23) 



{njaj)fo{d)+wfi{dy 

where all FPs which are responsible for the singularities are factorized. 



4 Sector decomposition of massless diagrams 

The sector decomposition of the FP integrals (13)-(16) consists of two steps. First a 
primary sector decomposition is performed that permits to eliminate the (5-function in 
such a way that the regions (20), which give rise to singularities, remain unchanged. In 
the second step, an iterated sector decomposition that permits to factorize all FPs that 
lead to singularities as discussed in the previous section is carried out. At the end of 
the sector decomposition those integrations that lead to 1/e poles and logarithms can 
be performed by means of standard formulas, which we derive in the appendices, using 
recursive subtractions of the type (21) and (22). 

In this section we describe the sector-decomposition technique to extract 1/e poles 
from massless^ diagrams [7]. A sector-decomposition algorithm to extract logarithmic 
singularities from massive diagrams is presented in Sect. 5. 



4.1 Primary sector decomposition 

As discussed in Sect. 3.3, our capability to isolate and simplify the singularities relies on 
the fact that they originate from regions of the type (20) . Therefore, when one eliminates 
the (5-function from the FP integrals (13)-(16) by performing the first integration, care 
must be taken that the structure of the singular regions remains unchanged.''' To this 
end, one can decompose the (/ — l)-dimcnsional hyperplane that is defined by the 6- 
function, into / primary sectors Vi, . . . ,Vf, defined as "Pj = {a | < for 1 < fc < / 
and J2j=i'^j = 1}) such that the sector Vi does not contain the hyperplane {a|aj = 0}. 
This yields 



G = J2Gi with Gi= dai [] da/t 1 - ^ a, ar[ 

i=l ■'^ ■'^ k=l \ j = l J 



(24) 



where we used the fact that 5 is a homogeneous function of degree — / in the integration 
variables, i.e. g{Xd) = g{a) as can be easily seen from (7). Now, each sector Vi can be 
mapped to the (I — l)-dimensional hypercube by means of the variable transformation 

with d^a = dai af^M^" V- (25) 





if l<k<i 


"fc = < 






if i<k<I 



''This technique is apphcablc also to massive diagrams, however, it permits to isolate only the 1/e poles 
and not the mass-singulax logarithms. 

^In particular, one cannot integrate one of the FPs, let's say a,, down to ai = 0. Otherwise, a singularity 
that might be located in the hyperplane {d\ai = 0} would be shifted to the hyperplane {a| '^j-^^ o-j = 1}. 
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Then, eliminating the 5-function through the ctj-integration one obtains 

"Ida, ^1 - / / 



(26) 



where 



gi{fi) = gy^^ = g{vi,---,Vi-iA,rii,---,m-i), (27) 

and using analogous definitions for Ui, J^i and Mi, 

gM = r(e + Le) ^."^'^"^^ , . (28) 

It is easy to see that this primary sector decomposition does not modify the structure 
of the singular regions (20). Indeed, each hyperplane {d\ak = and J2j=i (^j = 1}) which 
can give rise to singularities, has been divided into the sectors Vi with i k, and there 
it has been mapped to the hyperplanes {f]\r]i = ak/cti = 0} with I = k for i > k and 
I = k — 1 for i < k. Therefore, all singularities still originate from regions of the type 
{vlVji = ■ ■ ■ = Vjn = 0} with 1 < n < / — 1. 

4.2 Iterated decomposition of UV and mass singularities 

Now we perform an iterative decomposition of the primary-sector integrals (26) into sums 
of subsector integrals, where all FPs rji that give rise to singularities can be factorized. 
The UV singularities associated with the polynomial lAi arc decomposed as follows: 

1. If Wi(0) = 0, then we proceed with steps 2-4. Otherwise no further decomposition 
is performed. 

2. We choose a set TZ = {rjj-^ , . . . ,rij^} with a minimal number r of FPs such that 

Ui{fj) = if r]j, = ... = Vjr = 0. (29) 

Here, in order to simplify the notation, we assume that jk = k such that the poly- 
nomial Ui can be written as 

r 

Ui{f!) = Y.nkUik{f]). (30) 

k=l 

3. The integration range is decomposed into r subsectors Si,...,Sr with Sj = {f]\r]k < 
rjj for 1 < < r}, so that 

Gi = Y^Gij with Gi,= drjj II dm II dm9^{v). (31) 
j=i k=i Jo l=r+l 

Each sector Sj is remapped to the unit cube using the variable transformation 

ifr + l<k<I-l 

Cj^k otherwise 



e,- ifk=j with d'-'fj=d'-'ief'- (32) 
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As a consequence, in the sector Sj the variable can be factorized by rewriting the 
polynomial Ui as 



Ki{f])= ^jUij ii) with Uij ii) = Uij {f]) + J2 ikUik iv) ■ 



(33) 



k=l 



The resulting sector integrals read 
where 

g,j{0 = T{e + Le) 



, -> -] e+Ls 



f-{L+l)e 



(34) 



(35) 



with6 J^ijd) = J^i{fj) andUijii) =J\fM. 

4. For each subsector integral Gij, we restart the decomposition from step 1. 

The iterative application of steps 1-4 gives rise to a tree-like structure. At each iteration 
the subsectors are divided into new subsubsectors, which have to be labelled with new 
indices 



Gi — ^ Gij^ , 



G. 



n 



G 



(36) 



At each decomposition new FPs are factorized and the resulting sector integrals have the 
general form 

Gij,... = £ d^-if p(f^«.i...+-^n...) g,.^ ji), (37) 

where 



9m...i0 = r{e + Le)- 



e+Le 



(38) 



and the factorized FPs are written as p(^ih -+^'^ih -^ using the shorthands (87)-(88). 

The iteration of steps 1-4 stops when only sector integrals with Uij-^...{^) ^ remain. 
Then, the same iterative decomposition has to be repeated for the mass singularities that 
are associated with the polynomial J^ij^... until 



Wiji...(0) 7^ and J^iji...(0) 7^ 



(39) 



in all subsectors. 

The convergence and the efficiency of the sector-decomposition algorithm depend on 
the choice of the sets of FPs 71 = {rjj^., ■■■ ,r)j^} in step 2, which determines the subsectors 

®In some cases the FP can be factorized also in !Fij{$,) or A/ij(^). 
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Si,. . . ,Sr in step 3. In principle, the choice of minimal sets TZ, i.e. sets with a minimal 
number r of FPs, permits to reduce the number of new subsectors that are generated at 
each iteration and guarantees a certain efficiency. However, there are in general different 
minimal sets TZ that ftilfil (29), and it is not clear how to choose between them in order 
to minimize the total number of subsectors at the end of the iterated decomposition. 
Moreover, in some cases it is possible that owing to an unfortunate choice of the minimal 
sets the conditions (39) are never realized, i.e. the convergence of the algorithm, as it is 
formulated above, is not guaranteed. The approach that we have adopted to avoid this 
problem consists in a random choice of the minimal sets TZ. This strategy is in general not 
optimal but has the advantage to be non-deterministic, in the sense that every time that 
a certain integral is computed it is decomposed into different sums of subsector integrals, 
and this allows for checks on the final result. 

4.3 Counting of singularities 

In the following we focus on the properties of a generic sector integral. To this end we 
define 

G{T,e,g,:F) = [\'-'ip(i''+'A ^^L^e (40) 

with 

~9{0 = T{e + Le) ^^/g (41) 

and 

W(0) 7^ 0, :F{6) 7^ 0. (42) 

The sector integrals (37) can be written as Gij^,,, = G{Tij^,,,,e, gij^,,,, J^ij^,,,) with these 
conventions. The dependence of G on T, e, g and has been kept in explicit form for 
later convenience, whereas the dependence on all other quantities is implicitly understood. 
Note also that e can be different from e in (40)-(41). 

The FPs that give rise to singularities in the sector integrals (40)-(42) can be identified 
by simple power counting. In principle all FPs in 5 = {(,k\Tk < — 1} yield singularities 
unless these are not compensated by the behaviour of AA(^) at ^ ^ (J. Therefore, at the 
end of the sector decomposition, we decompose A^(^) into monomials with respect to the 
FPs in S and we treat the contribution from each monomial as a different subsector, such 
that the FPs S can be factorized and the resulting numerator J\f does not depend on them 
anymore. 

Now, those FPs that do and those that do not give rise to singularities can be immedi- 
ately identified. If we rename them as y = (yi, . . . , ym) and x = (xi, . . . ,xi), respectively, 
with I + m = I — 1, then the contribution from each subsector assumes the general form 

G ^ G{f,erg,:F) = /'d'f P(f^>-) fd^yPiyP^-'-') /^^'f^^^, (43) 

T{x\ y) 
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with 



and exponents 



g{x;y) = r{e + Le) 



A = {Ai,...,Ai), 
b = {bi,...,bm), 



M{x) 



U{x;y) 



(44) 



a = (ai, . . . ,a/j, 
/3 =(/?!,..., /3^). 



(45) 



By construction, all exponents (45) are integer numbers with Ai > —1 and bj > 0. The 
exponents A, b and a, {3 are related to T and r, respectively, in an obvious way. As men- 
tioned before, we have decomposed J\f{x, y) into monomials with respect to the FPs y and 
we treat the contribution from each monomial as a different subsector. The integrations 
over the m FPs y give rise to singularities up to the order l/e"*. 



4.4 Integration 

The integrations of the FPs yi, ■ ■ ■ ,ym for a generic integral of the type (43) are performed 
in Appendix B. This is done through recursive subtractions of the type (21), which permit 
to extract the leading and next-to-leading singularities, i.e. the poles of order l/e"* and 
l/e'^"^. The result for the integral (43) is obtained by applying (94) to the integrand 



ho{y) 



J d''xho{x;y) with ho{x;y) = p(^x 



A+de 



gix;y) 



^{x; y) 



e+Le ' 



Expanding = ho{x] y) in e as 



OO i 



^0 = E TT^o 



j=0 



and similarly'' T{e + s) and Af we obtain 



^ NLL 



VIA" f'.i 



d'f<'D|4°^(f;0) 



+ £ 



D|/iW(x;0)+^AA^'^/i^"^(f;0) 



i=l 



to NLL accuracy with 



'■0 

,(0) 



P(x" 



+ -^ + (L + l)lnU-LlnJ^. 



'^If e < 0, then the leading term of r(e + e) = X^jlo ^^^-'ie)£^/j\ contains a 1/e pole: 

r(°\e)=^^- for e<0. 
(-e)! e 



(46) 



(47) 



(48) 



(49) 



12 



The derivative operators and the subtraction operators A/ are defined in (95) and 
(96), respectively. These latter are associated to those next-to-leading contributions that 
remain after subtraction of the singularity from the yi integration. For instance, in the 
special case 6 = where all y integrations are logarithmically singular, we simply have 

A°'^4°) (f ; 0) = /' ^ [4°) (x ; 0, . . . , y,, . . . , 0) - (f ; 0)1 . (50) 
Jo Vi ^ J 

By construction, the remaining yi integration in the A^'* term is finite and can be per- 
formed either numerically or analytically. Also the Z-dimensional integral over the FPs x 
in (48) is convergent. The coefficients of the leading and next-to leading poles have thus 
been expressed in terms of convergent integrals with dimension I and I + 1. For negative 
kinematical invariants these integrals are well-defined. For positive kinematical invariants 
the analytic continuation is provided by the infinitesimal imaginary part e contained in J^. 



5 Sector decomposition of massive diagrams 

Let us now consider the sector decomposition of massive FP integrals (13)-(16) assuming 
a hierarchy of energy and mass scales as in (2). In order to extract the singularities that 
arise as combinations of 1/e poles and logarithms oi w = M^/s, which appear in the 
asymptotic limit w <C 1, we proceed as follows. 

First, wc perform the primary and iterated sector decompositions as described in 
Sects. 4.1-4.2. As a result, the sector integrals assume the form (40)-(42), where owing to 
(42) all FPs which lead to 1/e poles are factorized. Then, in each sector, we decompose 
the polynomial as 

HO=foii)+whii)-ie, (51) 

where /o(0 /i(^) are the contributions associated with the high-energy and mass 
parameters, and depend linearly on the corresponding ratios Sj/s and M^/M"^, respec- 
tively. If /o(0) = 0, there can be mass singularities that are regulated by the mass term 
wfi{S^) giving rise to logarithms of w in the asymptotic limit. In this case, we perform 
an additional iterated decomposition analogous to that of Sect. 4.2 until all FPs that lead 
to /o(0) = are factorized. At the end the original integral consists of a sum of sector 
integrals of the type 

m 



where 



-gii) = r(e + Le)—^^^, (53) 

and 

Hi) = Pii') /o(0 + wfiii) - ie. (54) 
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After the last sector decomposition we have 

W(0)7^0, :F{6)^0, and /o(0) 7^ 
in all subsectors. Without loss of generality we can assume that^ 



f = (ti, . . .tp,0, . . . ,0) with ti>0 for 



1, 



,P, 



(55) 



(56) 



i.e. that the FPs that have been factorizcd in (54) arc the first p FPs, where < p < / — 1. 

Simple power counting indicates that the set of FPs which can give rise to 1/e poles or 
mass-singular logarithms is given by 5 = {^k\Tk < —1 or — etk < —1}. As discussed in 
Sect. 4.3, the numerator A^(^) has to be decomposed into monomials with respect to the 
FPs in S, and each monomial has to be treated as a different subsector where the FPs in 
S can be factorized. 



5.1 Transformation of non- logarithmic singularities 

Let us discuss the behaviour of the integrations over the FPs ^1 , . . . , in the asymptotic 
limit w ^ 0. If Tfc > — 1 for 1 < A; < p, the ^jfc-integrations are convergent in D = A 
dimensions. For w they behave as 



dCk 



^k +w 



^^tkK^k 



tk\e—di. — l 



^k +W 



1 if 4 < 
ln(u;) if 4 = 
w-'^'^ if 4 > 



(57) 



where we have introduced the degree of singularity 

~ Tk + 1 



dk 



tk 



(58) 



for k = l,...,p. Note that if 4 > 0, the ^fc-integral yields a non-logarithmic mass 
singularity of order w"'^'' for > —1. Similarly, for < — 1 one obtains singular 
contributions of order w~'^''e~^ in D = 4 — 2e dimensions. Independently of the values 
of T^., these non-logarithmically divergent integrals can be easily expressed as a linear 
combination of logarithmically-divergent integrals with dk < for k = 1, . . . ,p. To this 
end, we define the maximal degree of singularity 



d := max (41 
i<fe<p 

and if d > we multiply the integrand in (52) by 

d 

= (wfi - ie) J2 



(59) 



wfi 



le 



r=0 



■d—r 



-/oP(f) 



(60) 



This yields 



G{T,e,g,:F) 



w 



r-=0 



{-lYG{Tr,er,gr,J^), 



(61) 



This requires a permutation of the FPs that is, in general, different for different subsectors. 
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where 



Tr = T + rt, er = e-d + r, gr = gifoYifi - ^e^^ ■ (62) 

As a result of the above transformation, the term w~'^ involving inverse powers of masses 
is factorized and the degree of divergence dk is reduced by d in the resulting integrals 
G{Tr,^r,gr,^)- This can be easily seen from 

[dr) ^ := - =dk-d. (63) 

Thus, the integrations with maximal degree of singularity dk = d become logarithmically 
divergent, whereas all other integrations with dk < d become non-divergent. 

Since the integrals resulting from (61)-(62) have the same structure as the original 
ones, we need to compute only integrals of the type (52)-(56) with 

dk<d<0 for k = l,...,p (64) 

in the following. Integrals with Tk < —1 for some 1 < k < p, where the corresponding 
^fc-integrations give rise to additional 1/e poles, do not need to be computed explicitly^ 
since, as we show in Appendix D, this kind of integrals can be eliminated by means of 
integration-by-parts identities. Thus, we can restrict ourselves to integrals with 

Tfc > -1 for k = l,...,p. (65) 

We note that, as a consequence of (64)-(65), 

d = ^ e > 0, (66) 

since d = implies that there is a A; with I < k <p and dj- = 0. Thus e = (T^ + l)/tk > 
owing to (58) and (65) . 



5.2 Classification and counting of singularities 

In the following, we consider sector integrals of the type (52)-(56) with (64)-(65). In 
order to extract the various 1/e poles and ln{w) contributions, we first classify the FP 
integrations in (52)-(56) according to their singular behaviour. 

• The FPs {^fe|l < k < p and dk = 0}, for which we can assume tfc > and Tk > — 1, 
are renamed as {zi, . . . , Zn}- As we see in Sect. 5.3, the corresponding n integrations 
give rise to singular contributions up to the order ln^{w). 

• The FPs {^j\p + l<J<-f — 1 and Tj < —1}, for which we can assume tj = 0, are 
renamed as {yi, . . . ,ym}- The corresponding m integrations give rise to poles up to 
the order l/e"*. 

®In practical calculations, wc never encountered contributions of this type at the NLL level. However, 
in general we cannot exclude their presence, especially in view of an extension of the algorithm to the 
NNLL level or beyond. 
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• The remaining FPs, {^i\Ti > —1 and if 1 < i < p. di < 0}, are called {xi, . . . ,xi}. 
The I integrations over these FPs are free of singularities. 

In general we have < l,m,n < I — 1 and / + m + n = /— 1. As a result of the above 
classification, the contribution from each subsector assumes the general form 



G = G{f,e,g,:F) = 



i+^e\ £l{x;y;z) 



J 



T{x;y;z) 



e+Le 



with 



g{x;y;z) = r(e + Le) 



AA(x) 



U{x;y;z] 



f-{L+i)e 



(67) 



(68) 



and 



r{x;y;z) = P(^f") ^if'') fo{x;y; z) + wfi{x;y; z) - ie, 



(69) 



where U{0: 0; 0) / 0, ^(6; 0; 0) ^ 0, and /o(0; 0; 0) ^ 0. The exponents that are associated 
to the FPs X and y, 



A = {Ai,...,Ai), a = {ai,...,ai), a = {ai, . . . ,ai), 

b = {bi,...,bm), 0={Pi,...,Prn), 



are integer numbers and satisfy 

Ai>-1, Ai>aie-1, 



ai > 0, 



bj > 0, 



(70) 



(71) 



by construction. The number n of FPs z depends on the maximal degree of singularity d. 
If d < then n = 0, otherwise n > and for the exponents of the FPs z, 



ce = (cie, ...,Cne), 7 = (71, . . . ,7„), 
which are integer numbers, we have 

Cfe > 0, Cjte > 0. 



(72) 



(73) 



Again, the exponents A, a, b, c and a, /3, 7 are related to T, t and r in an obvious way, 
M{x, y, z) has been decomposed into monomials with respect to the FPs y and z and the 
contribution from each monomial has been treated as a different subsector. Note that 
n > 0, or equivalently d = 0, implies e > 0, as already observed in (66). 
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5.3 Integration 

Before we proceed, we rewrite (67)-(69) as 



11 1 P( ^'^^■"-'-+7^ ) 

G = [ d'x f iTyPiyP'-^'-^) f z ho{S- y; z) --—^ ^ 

^0 ^0 ^ ''Jo \P z'^ +wh^ x:v:; 



with 



and 



[P{z<^+wh,ix;y;z)Y+''' 
9{x;y;z} 



hi{x;y;z) 



foix;y;z) - ie 

/i(,r: lj:z) - u 



e+Le 



fo{x;y;z)- ie P(x«) 



(74) 



(75) 



(76) 



The singular integrations over the FPs y and z for a generic integral of the type (74) are 
performed in Appendix C to NLL accuracy. The resulting 1/e poles and ln('i/;) contribu- 
tions can be read off from (113). Expanding ho = ho{x; y;z) in e as 



oo j 



(77) 



and similarly r(e + e) and N, we obtain 



G 



NLL 



P(p)P{c. 



(-irE 7 



p>0 



d'f i Fp{f)j£hf{x; 0; 0)L„+p {w) 



+ 



e{n - l)Fp(f)D| [4°) (f ; 0; 0) [in (hi{x- 0; 0)) + Cg 



-E^.-^p(%i)4'd|4'^(^;0;0) 



where L^j is defined in (83) and 

C = rW(e)p(x 

_ rW(e) 



d| (x; 0; 0) + ^ ^3^l^fhf (x; 0; 0) 



(78) 



tA—ea 



Ar{o) 



(/o-ie)'z://' 



^fmE + 1" (^(-"-"')) + ^ + (i^ + 1) InW - Lin (/o - ie) . (79) 



;,(0) r(0)(e) 



The derivative operator and the subtraction operators A^' and AV are defined in (95) 
and (96) and (108), respectively. The constant Cg is the sum defined in (85). The vector 
r = (ri, . . . , r„) has components 



7i 



-L, 



(80) 



17 



and ry] and Fp{r) arc defined in (89) and (90), respectively. 

As observed at the end of Appendix C, the result (48) for massless integrals, can be 
obtained from the general result (78) for massive integrals in the special case n = 0. In 
the massive case, for the integrals (67)-(69), where the integrations over the FPs z are 
logarithmically singular, the m + n integrations over the FPs y and z give rise to singular- 
ities up to the order e~'^^^ \n^~^^ {w) with p > 0, as one can see in (78). Instead, for the 
asymptotic behaviour of non-logarithmically singular integrals we observe the following. 
For an integral (52)-(56) with d > and (65), which is expressed in terms of logarithmi- 
cally singular integrals of the type (67)-(69) using (61)-(62), each term in the sum (61) 
yields a contribution of order e~"^^^w~'^\n'^^^{w) with p > 0. However, from (62), (75) 
and (78) one can see that these leading contributions are independent of r and thus cancel 
in the sum (61), i.e. to LL accuracy 

r=Q ^''^ 

w-^Gif, ~e - d, fr^-g, ^) E f"^) i-iy = 0. (81) 

r=0 ^''^ 

Consequently, only the subleading contributions proportional to Cg^ survive, i.e. contri- 
butions of order 

^-m+pyj-di^n+p-lf^^j with p>0. 

The result (78) is expressed as Z-dimensional integral over the FPs x. Analogously 
to the massless case, the next-to-leading contributions involve one additional integration 
over the FPs yi or zj, which appears in the subtraction operator A^'* or A^'-' , respectively. 
Thus, the coefficients of the LLs and NLLs have been expressed in terms of convergent 
integrals with dimension I and I + 1. For negative kinematical invariants these integrals 
are well-defined. For positive kinematical invariants the analytic continuation is provided 
by the infinitesimal imaginary part ie contained in hg and hi . 



6 Discussion 

The sector-decomposition algorithm that we have presented in Sects. 4 and 5 permits to 

reduce arbitrary massive (or massless) FP integrals to a sum of sector integrals of the type 
(67)-(69), where all FPs that can lead to UV or mass singularities are factorized. Each 
sector integral gives rise to a tower 

J oo 

E E a,,ke'ln^^Hw), (82) 
j=0 k=-j 

involving 1/e poles and logarithms. As discussed in Sect. 5.2, the FPs yi,---,ym and 
zi,...,Zn, which give rise to these singular contributions arc easily identified by power 
counting. The maximal total power J of the resulting poles and logarithms is given by 
the number of these FPs y and z, i.e. J = m + n. This simple counting permits to select 
the sector integrals that give rise to the LLs and NLLs by requiring that J = 2L and 
J > 2L — 1, respectively, at L-loop level. Here wc have assTuncd that the overall factor 
r(e-|-L£) in (15) does not yield 1/e poles since the LLs and NLLs originate from diagrams 
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involving mass singularities, i.e. diagrams with an exponent e > for the mass-dependent 
denominator .F in (15). 

The singular contributions have been extracted from the integrals over the FPs y and 
z in (67)-(69) using the subtraction technique discussed in Sect. 3.3 and Appendix C [see 
(109)-(110)]. Each of these J integrations has been split into a singular part that has been 
solved analytically (giving rise to a pole or a logarithm) and a non-singular part that has 
not been integrated in (78). The contributions with total power j in the resulting tower 
(82) originate from the combinations of / singular parts and J — / non-singular parts 
with j < j' < J. Their coefficients are thus expressed as integrals over (J — j') remaining 
FPs Ui or Zk- Finally, together with these FPs, also the FPs xi, . . . ,xi, which do not give 
rise to singularities, remain to be integrated. In the general result (78), the coefficients of 
the leading (j = J) and next-to-leading {j = J — 1) contributions for the generic sector 
integral (67)-(69) are expressed as (/ — 1 — J)- and (/ — J)-dimensional integrals, where 
/ = l + m + n + 1 is the number of internal lines of the Feynman diagram. These remaining 
integrations are free of UV and mass singularities by construction. 

For the L-loop corrections to processes with E external legs, which involve a maximal 
number I = E + 3L — 3 of internal lines, the coefficients of the LLs (j = 2L) and NLLs 
{j = 2L — 1) consist of (E + L — A)- and (£' -|- L — 3)-dimensional integrals, respectively. For 
the most complicated Feynman diagrams that we have computed up to now [18,19], which 
correspond to 3- and 4-point functions at two loops, the NLLs require 2- and 3-dimensional 
integrals, respectively, and such integrals turn out to be sufficiently simple to be solved 
analytically by standard computer algebra programs such as Mathematica[27]. 

To our knowledge, the algorithm we have presented is the only existing tool that 
permits to extract NLLs from arbitrary Feynman diagrams in a completely automated 
way. At present, the accuracy is limited only by the fact that the divergent integrals have 
been solved in NLL approximation. In order to extend the algorithm beyond the NLL 
level it is sufficient to solve these divergent integrations, which correspond to a well-defined 
class of standard integrals resulting from the sector decomposition, with higher accuracy. 

The algorithm has been checked against all one-loop tensor integrals given in Ref. [24] 
and the two- loop master integrals of Ref. [5] . Its first applications were the two-loop cal- 
culations of the angular-dependent subset of the electroweak NLLs for arbitrary processes 
[18] and the complete set of NLLs for the electroweak-singlet massless fermionic form fac- 
tor [19]. In the future, it can be applied to derive the two- loop electroweak logarithms 
for more complicated processes in order to check the existing resummation prescriptions, 
which still rely on arguments based on symmetric gauge theories. Alternatively, it can be 
used to extend the approach that has been used in Ref. [28] from the one- to the two-loop 
level, in order to derive the two-loop NLLs for generic processes within the spontaneously 
broken electroweak theory. 

7 Conclusions 

At energies far above the electroweak scale, ^/s S> Mw, the electroweak radiative cor- 
rections are dominated by large logarithms of s/M^^. These logarithms originate from 

^"The additional integration for the coefficients of the NLLs is Eissociated with the subtraction operators 
A^'" and A°J' defined (96). 



19 



ultraviolet and mass singularities. In this paper, we have presented an algorithm that per- 
mits to extract these singularities and the corresponding large logarithms from arbitrary 
multi-loop Feynman integrals by using sector decomposition. We have elaborated the algo- 
rithm for the case of a single mass scale and a single energy scale but allowing for various 
different mass and energy parameters. This permits, in particular, to compute higher- 
order next-to-leading logarithmic electroweak corrections for processes involving various 
kinematical invariants of the order of hundreds of GeV and masses Mw ~ Mz ~ Mh ~ mt 
of the order of the electroweak scale in the approximation where the masses of the light 
fermions are neglected. 

We have provided explicit formulas for the extraction of the leading and next-to- 
leading mass singularities. At the next-to-leading level, the algorithm has been success- 
fully checked against one-loop and two-loop results available in the literature, and first 
applications have already been published. This method will be useful for the calculation of 
next-to-leading logarithmic electroweak two-loop corrections and for the check of existing 
resummation prescriptions at this level. 

The method can be extended beyond the next-to-leading level. To this end, a well- 
defined set of mass-singular integrals has to be solved in the required approximation. 
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Appendix 

A Notation and conventions 

For logarithms we introduce the notation 



Ln{w) := 



e{n) 



In' 



(83) 



n! 



The 9 and the 5 functions with integer arguments are defined as 




for n < 



for n > 




for n = 



for n 7^ 



(84) 



In order to keep our results as compact as possible it is useful to define 




(85) 
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Note that Ce = *(e) -|-7e, where *(e) = r'(e)/r(e). 
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Here and in the following all sums have to be understood as 



E=^(ii-io)E- (86) 



J=JO 3=30 



Similarly, we use the convention IljLjo • • • = 1 Jo > Ji • 
For vectors z = {zi, . . . , Zn) we introduce the shorthands 



rl rl n n 

/ dV:= / IJd^i' P{z):=^z^. 
Jo Jo fj, fL\ 



(87) 



In the special case dim(z) = n = the above expressions have to be understood as 
Jq dPz = Jq Y[^=i dzj = 1 and P{z ) = 11^=1 = 1- Powers of vectors z = {zi, . . . , Zn) 
with vector exponents a = (ai, . . . , an) or scalar exponents b are defined as 



The (n — 1) -dimensional vector resulting from the omission of the jth component of an 
n-dimensional vector z = (zi, . . . , Zn) is denoted as 



\3\ 



(Zl, . . .,Zj-i,Zj+i, . . .,Zn). (89) 



For vectors r = (ri, . . . , r„) with real components and dimension ra > we introduce the 
function 

p / n \ n 

F,{r):=eip) Yl Hp-J2li)U^f (90) 

qi,...,q„=0 ^ j=l ^ j=l 

involving all monomials of degree p. In the special case n = we have Fp{r) = 5(j)). Note 
also that -Fo(r) = 1 for n > 0. 



B Massless integrals 

In this section we compute FP integrals of the type (43), which result from the sector de- 
composition of massless loop integrals. In particular, we perform the singular integrations 
over the FPs y to NLL accuracy, i.e. including the leading and next-to-leading 1/e poles. 
Let us first consider one-dimensional integrals of the type 



['dyyP'-'-'ho{y) (91) 
Jo 



with P real and 6 > integer. We assume that the singular behaviour of the integrand at 

y ^ has been isolated in the term yl^'^~'^~^ such that /iq is an analytic function free of 
poles and finite at y = 0, i.e. /io(0) 7^ 0. The 1/e pole, which originates from an end-point 
singularity at ?/ = 0, can be isolated by adding and subtracting from /io(2/) the term 



/io(2/) = E/D^/io(0) with D^/io(yo,^- ^ 



k=o ^- 



(92) 

y=yo 
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To NLL accuracy we obtain 

fi 



/ = / dyy 
Jo 

b 



E 

1 



1 



ho{y) + / dyy 





,/3£-6-l 



ho{y) - hoiy) 



b-l 



k=0 



y" 



ho{y)-hoiy) +0{e). (93) 



As usual in dimensional regularization, divergent integrals with b> are defined via ana- 
lytic continuation. The generalization to m-dimensional integrals is obtained by recursive 
application of (93). This yields 



d"'9P(?*-*-') hoiy) 



ill- 1 _ - /I 



m— 2 



(94) 



where the derivative operator and the subtraction operator A^'* are defined as 



^ho{yo) = 



n 



hi 



Qbl Qbr, 



■hoiy) 



(95) 



y=yo 



and 



Ay/io(yo)= E 

k 

dyi 



y ^ h- — h- 



^hoiyo) 



bi=ki 



+ / ~b^\ I'^ghoiyoi, . . . ,yoi-i,yi,yoi+i, ■ ■ ■ ,yom) 
Jo y-'^ [L 



bi=0 



ki=0 



b-l — k-i 



(96) 



The m integrations over the FPs y in (94) yield leading and next-to-leading singularities 
of order and l/e"*"^, respectively. Note that the subtracted terms (96), which enter 
the coefficients of the next-to-leading poles, involve one-dimensional integrals over one of 
the FPs yi. These integrals are free of the singularity at = since the corresponding 
terms within the curly brackets in (96) are of order y^*^^. 



C Massive integrals 

In this section we compute FP integrals of the type (74), which result from the sector de- 
composition of massive loop integrals. In particular, we focus on the singular integrations 
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over the FPs y and z. To be specific, we consider the generic integral 

/ d^yPiy^^-'-') d^zh,{y-z) ^ . ^m^+l. - (97) 

Jo ^ 'h {P{z^ ^whx{y\z)\ 

We assume m, n and e integer with m, n > and e > if n > 0.^^ The exponents 
(5 = {Pi,. . . , Pm) and b = (6i, . . . , bm) are integer numbers. Moreover, we assume 6j > 
such that the integrations over the FPs y = {yi, . . .ym) give rise to a singularity of order 
1/e™. The exponents 7 = (71, . . . ,7^) and ce = (cie, . . . ,Cne) are integer numbers with 
Cj > 0. The integrations over the FPs z = {zi, . . . Zn) give rise to mass singularities that 
are regulated by the mass-to-energy ratio w = M^/s. 

All mass and UV singularities appear as end-point singularities at y — or z — 0. 
As a result of sector decomposition, all FPs that are responsible for such singularities 
are isolated in the terms p(jjf^^~^~^^ and P^z^^. Thus, /iq and hi, which arc analytic 
functions corresponding to (75) and (76), are free of end-point singularities and finite at 
y — 0, £ — > 0, i.e. /io(0, 0) / and hi{0,0) 7^ 0. We also assume that possible poles of 
ho and hi that approach the integration contour as ie — ^ can be avoided by appropriate 
contour deformations. 

The integration is performed in the asymptotic limit 

< u; < 1 (98) 

to NLL accuracy, as explained in Sect. 2. In the following we present some intermediate 
results which correspond to special cases of (97) and have been used to compute the above 
integral in the most general case: 

• Let us start with the special case ho = hi = 1 and m = 0. Here we define 

where we include combinations of logarithms which originate from the e-expansion 
of the integrand in (97). For q = {qi, . . . , qn) with integer qj > we obtain 

n 

Bl^w)""^^ {-lT[hf,{w) + C,hf,_i{w)\ with h = n + Y,(lj. (100) 

where Cg is defined in (85). For n = the result (100) is trivial, since 

Bl{w) = j^^^=l + 0{w). (101) 
For n > 1 and e = 1 we have a polylogarithm 

<,» = (-l)"-'^+^Li,(-l) 

= (-1)" hf, {W) + O [l.n-2 {W)\ +0{W). (102) 



^^As noted at the end of Sect. 5.2, for n > we need only the case e > 0. 
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For n> 1 and e > 2, (100) can be proven by induction using the recursive relation 



^n,(,i,...,5„)(«^) = ^;(,,,...,5„)(«^) + ^ 



B 



n,{qi,...,q„-i,q„-l)^ ' 



which follows from partial integration in Zn- 
• Using (100) one can easily prove that 



(103) 



(Piz')) (-l)"i^p(f) [L„+p (^) + CeL„+p_i («;)] 



(104) 



for any real vector r = (ri, . . . , r„) and integer number p, where -P'p(r) is the combi- 
natorial factor defined in (90). 

In the special case /iq = /ii = 1 and 6 = we obtain 



/'d'».-p(?*-')rdv- 



NLL 



P(/3)P(C)] (-l)»5]F/-' 



p>0 



[Ln+p (w^) + C'eLn+p-l (w^)] 



where the components of the vector r = (ri, . . . , r„) read 



L. 



(105) 



(106) 



The integration of the variables yi, . . . , in (105) is trivial and gives rise to a pole 
of order 1/e"*. The remaining integrand has to be expanded up to the needed order 
in e. Then, after neglecting irrelevant terms of order u;, it can be integrated with 
the help of (104). 

In the special case m = 0, including logarithms as in (99), we obtain 

P{z~^-^) 



I'd^'zhoi 
Jo 



P{z) + whiiz) 



n ^^{Zj) (-1)" /io(0)L^i {wh.iO 



+ 



/io(0)Ce-^%)A°'^/io(0) 



Lfi-i (w) 



(107) 



where n = n + J2j=i Qj analogously to (96) the subtraction operator A^'-' is 
defined as 



A^'^/io(^o) = / —[ho{zoi,...,zoj-i,Zj,zoj+i,...,zon)-ho{zo)]. (108) 



1 dz. 



Zj 
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The above result has been obtained by means of the subtraction technique discussed 
in Sect. 3.3. Each integration has been spht into singular and non-singular parts as 



/ dzjI{ho,hi;z) = [ dzjI{ho,hi; z) 
Jo Jo 

+ J dzj[l{ho,hi;z)-I{ho,hi;z)\f^.^f^ 



hi=hi{zi,...,Zj-l,0,Zj+l,...,Zn) 



{zi,...,Zj-l,0,Zj+l,...,Zn) 



, (109) 



where /(/lo, hi;z) represents the integrand on the left-hand side of (107). The recur- 
sive application of this subtraction to all n integrations gives rise to various terms, 
which consist of combinations of I singular parts and n — I non-singular parts with 
< I < n. The LL and NLL contributions correspond to the terms with I = n,n—l, 
and read 



hi=hi(0) 



C d''zI{h^M;^''= f^d''zI{ho,h,; 
Jo Jo 

+ E / dV[7(/io,^i;^U,=h,(o,...,o,^,-,o,...,o) - ^(^o,^i;^U.=h,(o)J • (110) 
i=i •'^ 

The first term on the right-hand-side of (107) originates from the first term on the 
right-hand-side of (110) and has to be expanded to NLL accuracy as 



Ln [whi{Q)) Lf, (w) + In (hi{Q))Ln-i (w) . 



(Ill) 



The second term on the right-hand-side of (110) gives rise to the A contributions 
in (107). Here the function S{qj) indicates that only terms with qj = yield a NLL 
contribution. 

• Using (107) it is easy to prove that for any real vector r = (ri, . . . , r„) and integer 
number p, 

£d^zho{^ ^^^'"'^ g Lp (^(-^"O) (-l)"|^o(0)Fp(r-)L„+p (wh,{6)) 



+ 



P{z)+whi{z) 



-'n+p—l 



(112) 



where f^j and Fp are defined in (89) and (90), respectively. 
Finally, combining (94) and (112) we arrive at the result 



■/3e-b-l 



d'^zho{y;z) 



[P{z'^ +whi{y;z^ 



1 e+Ls 



NLL 



p(/3) Pic)]'' i-ir E (7)"'' i^p(r-)D|4°)(0;0)L„+^ (w) 



p>0 
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+ 



e{n - l)Fp(r)D| 0) [in (/ii(0; 0)) + C, 



-E9^p(%])aJ'^d|4°)(0;0) 



d|/.W(0;0)+^AAJ/4°\0;0) 



-'n+p 



-iH 



where 



(i) 



j=0 



is the e-cxpansiori"'^^ of the function ho = hQ{y;z). The vector f 



[n . . . ,r„j, ry], 



(113) 



(114) 



and Fr, 



are defined in (106), (89), and (90), respectively. Moreover, the derivative operator and 

the subtraction operators A^'* and A^'"* are defined in (95), (96) and (108), respectively. 
In particular, we have 



A°>D^-4°)(0;0) 







DC4°)(0;0,...,z,-,...,0)-D|/i^"^(0;0) 



b j,(o)('n. n^ 



(115) 



The m + n integrations over the FPs y and z in (113) yield leading and next-to-leading 
singularities of order £-"^+p ln"'''^(u') and Itl"'~^^~^{w), respectively. Similarly to the 

massless case, the subtracted A terms which enter the coefficients of the next-to-leading 
poles, involve convergent one-dimensional integrals over one of the FPs yi or Zj. 
We note that in the trivial case n = the integral (97) simplifies to 



hojy) 
[1 + whi '^-^^~^+^' 



= £d"'yP(y 



ho{y) + 0{w). (116) 



Since no singularity is regulated by the mass term w, in the asymptotic limit u; — > (116) 
corresponds to the massless integral that is obtained by setting w = Q m. the integrand. 
Indeed, in the special case n = 0, the general result (113) for massive integrals is equivalent 
to the result (94) for massless integrals, as one can easily verify using Fp{f) = 5{p) for 
dim (r) = n = 0. 



D Factorization of poles from massive integrals 

In this section, we consider again the generic sector integrals (52)-(56), which have been 
computed in Sect. 5 assuming (65), such that all singularities resulting from the integra- 
tions over the FPs ^i, • • • , Cp, i.e. the FPs that have been factorized in (54), are regulated 
by the mass term w in (54). Here we discuss the more general case where < — 1 for 
some 1 < k < p and the corresponding ^^-integrations yield additional 1/e poles. As we 
have already anticipated, this kind of integrals can be eliminated by means of (recursive) 
integration by parts in the variables ^fc, which permits to increase the exponents Tfe < — 1 
until the 1/e poles are factorized and only integrals respecting (65) remain. 

^^The leading term of this s expansion can be of order 1/e*. 
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In the following we consider sector integrals (52) (56) with (64) assuming that possible 
non-logarithmic singularities have been eliminated as described in Sect. 5.1. An integration 
by parts in the variable yields 

where the ^-dependence of g and T is implicitly understood. The 5(1 — ^k) contribution 
on the right-hand side represents the boundary term at = 1, whereas the boundary 
term at = vanishes. -"^^ The remaining ikd^,. contribution being proportional to has 
the desired property to increase the power Tj. associated with the FP ^j.. However, as one 
can see from 

the degrees of singularity dj = e — {Tj + l)/tj of the various FP integrations 1 < j < p can 
increase for the term proportional to ttJ^fc3^^/i, where the increase of e is not compensated 

by the factor in tiie numerator. In order to preserve (64), this must be avoided by 

rewriting 

= _ pf^^ /J , (119) 

fi-ie ^ \ / J 

in (118). The integration by parts in the FP combined with (119) results in a sum of 
three contributions, 

G{T,e,g,J^) = Tfr—rr Gi'^r,er,~9r,:F), (120) 

that have the same structure as the original integral with 

ei = e, f 1 = f , gi = 6{1 - ^k) 9, 
62 = 6+1, T2 = f + t, g2 = ie + Le)g 



/i -le 



63 = e, (n),=Tj + 5jk, h = {e + Le)S^~g-d^J. (121) 
^ /i-ie 

The exponent associated with the FP grows in all terms, i.e. (T^)a; > Tfe, apart from 
the boundary term G{Ti,ei,gi,J^), where the ^fe integration is irrelevant. Thus, in order 
to eliminate all integrals that do not satisfy (65), it is in principle sufficient to iterate the 
integration by parts to all FPs with < —1. However, since our aim is to reduce 
by one the number of integrations every time that we extract a singularity, care must be 
taken when 1/e poles are extracted through the overall factor l/{Tk + l + Te) in (120), i.e. 
every time that we perform an integration by parts in a FP with = —1. In practice, 

^"'This holds for arbitrary values of Tk within dimensional regularization, i.e. also for Tfc < — 1 where the 
boundary term at ^fc = is divergent for e = 0. 
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as we show below, it is possible to restrict the use of (120)-(121) to the cases where e = 
and d^f^g = if T/t = —1, such that G{Tr,er,gr,^) = 0{e) for r = 2,3 and the only 1/e 
singularity results from the boundary term G{Ti,ei,gi,J^). 

Let us now outline our general strategy to eliminate all contributions that do not 
satisfy (64) and/or (65). Starting from a generic sector integral (52)-(56) we proceed as 
follows: 

1. If non-logarithmic singularities are present, i.e. d > 0, we eliminate them by means 

of (61)-(62), such that the resulting integrals satisfy (64). 

2. If there is a FP with 1 < k < p and < —1, this FP is integrated by parts using 
(120)-(121). This step is repeated until all integrals satisfy > —1 ioi 1 < k < p. 
Here no 1/e pole appears. 

3. For the remaining integrals involving a with 1 < k < p and = — 1 we have 
e = dk<d<Oasa result of step 1 and owing to (58). If e < 0, then 

-p-{e+Le) ^ p-Le | ^p((^ /q] + 0(«;)| , (122) 

and the contributions of order w can be omitted. Thus, we can rewrite the corre- 
sponding integral as 

G{f, e, ~g, T) = G{f - et, 0, gf^', T), (123) 

and the integral on the right-hand side, with exponents T ^ T — et, already satisfies 
(65) as a consequence of step 2 {Tk > — 1) and e < 0. 

4. Finally, we have to consider the integrals with e = and a with 1 < k < p and 
Tk = —1. In this case, we split g as 

m = hii) + ^9kii) with gkiO = 5(6, • • • , ^k-iA Cik+i, • • • , C/-i) (124) 
and treat these two contributions separately, i.e. 

G{f,e,~g,^) = G{f,e,gk,^) + G{f,e,A~gk,^). (125) 

On the one hand, for the contribution G{T, e,gk,J^) we can safely perform an inte- 
gration by parts in since, as discussed above, owing to d^^^gkiCj = and e = 0, 
only the boundary term effectively contributes to the resulting 1/e singularity. On 
the other hand, the remaining contribution is free of singularities owing to 

lim 4-iA5,(a = %5(aL 0' ^12^) 

and can thus be written as 

G(f, e, A^fc, ^) = G{f', e, ^^'^h, ■^), (127) 

where T' = (Ti, . . . , Tjfc_i, 0, T^+i, . . . , Tj^i). The power of has been reduced 
in both contributions without generating unwanted 1/e poles. At the end, if other 
FPs with 1 < j < p and Tj = — 1 are still present this step has to be iterated until 
all integrals fulfil (65). 
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